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Abstract:  A  user-friendly  ‘ divide-and-conquer'  algorithm  is  presented  for  finding  all  the  self- 
intersection  points  of  a  parametric  curve  in  the  Bernstein- Bezier  representation.  The  underlying 
idea  of  the  algorithm  is  to  deal  with  the  Bezier  polygon  instead  of  the  curve  description  itself. 
By  alternately  subdividing  the  Bezier  polygon  and  estimating  the  self-intersection  regions  the 
self-intersection  points  are  finally  approximated  by  straight  line  intersections  of  the  refined  Bezier 
polygons.  The  algorithm  also  calculates  the  parameter  values  of  the  self-intersection  points.  In 
addition  to  the  convex  hull  and  the  approximation  property  of  the  Bezier  polygon  the  working 
of  the  algorithm  is  based  on  a  very  intuitive  angle  criterion. 


0.  Introduction 

For  two  explicit  given  curves  f^(x)  and  intersection  points  of  /[(ar)  and  fi(x)  can  be  calcu¬ 
lated  Using  numerical  methods  like  Newton's  method  by  rewritting  the  problem  as  that  of  fmding 
the  roots  (zeros)  of  the  function  F(,x)  =/I{.’c)  -  fjix).  If  the  equation  of  one  curve  is  given  in  implicit 
resp.  e.xpltcit  form  and  the  other  in  parametric  form,  we  can  substitute  the  parametric  form  into  the 
implicit  resp,  explicit  equation.  The  (usually  non-linear)  equation  we  obtain  can  be  solved  by- 
Newton  s  method  again.  If  both  curves  are  given  implicitly  as  (non-linear)  functions  f\[x,y)  and 
i  of  .r  and  y  or  as  parameterized  curves  .e,  =  y,  =>.(0  yj  =  y2(T)  we 

have  to  solve  the  two  equations  /(jr. ;,  )  =  0  and  /^(.x.y)  =  0  resp.  .rdO  -  .X;(t)  =  0  and 
}  i!)  -  y-jr)  =  0  simultaneously,  what  can  also  be  done  by  Newton's  method  [faivc.Prait  'fJJ. 
geometncaly  based  modification  of  the  methods  described  by  Faux  and  Pratt  to  calculate  the 
intersection  points  of  two  parameterized  curves  was  given  by  Hoschek  in  [Hoschek  -Pi].  Hoschek  s 
method  works  also  for  the  problem  of  calculaimg  the  self- intersections  of  a  curve.  Self-intersections 
of  a  curve  con  appear  for  e.xample  as  boundanes  of  loops  of  parallel  curves,  often  called  olTset 
curves.  \_.Arnold  86],  [Farouki  '•Jd],  [Hoschek  -Si.  '87],  [Klass  '53],  [Lyche,.\iti>rken  -5^. 

[Tiller. Hansert  8T].  For  the  loop  removal  the  self-intersection  points  have  to  be  detected  For  ra¬ 
tional  curves  this  can  also  be  done  by  algebraic  methods  which  have  been  mtroduced  in  the  area 
of  CAGD  by  Sederberg,  Goldmann  and  .-Vnderson  They  desenbed  in  [Sederberg  ’84],  [Sederberg 
el  al.  "id,  85]  and  [Goldmann  '5i]  a  method  ot  classical  algebraic  geometry  for  solving  the  curv  e- 
curve  mtersection  problem  for  rational  planar  and  non-planar  curves  and  in  [Sederberg  et  al.  .Si] 
a  method  for  finding  the  double  points  and  by  this  the  seil'-intersection  points  of  planar  rational 
cubics  (see  also  [Salmon  18~9].  [Htlion  Ji],  [If  u.'/itT  it)]). 


n 


In  CAGD  the  B-spline- Bezier  representation  of  curves  is  very  popular  and  therefore  it  is  of  im¬ 
portance  to  have  (self-iintersection  algonthms  for  this  type  of  curve  representation  too.  so  that  no 
conversion  of  the  curve  description  iDannenberg.Swacki  '55],  [Hoschek  'S~]  is  necessar.. 
Curve-curve  intersection  algonthms  for  B-spline-Bezier  representations  have  been  desenbed  by 
ILane  ei  ai  50],  [_Cohen  el  al.  50]  and  for  quadratics  by  [Yang  et  al.  'Sb],  '^'ang  calculates  the 
mtersection  points  using  an  algebraic  method  while  the  algorithms  of  Lane  and  Cohen  are  subdi¬ 
vision  algonthms  taking  mto  account  the  geometnc  relationship  between  the  curve  and  its  defining 
control  polygon.  Pure  subdivision  algonthms  are  very  time-consuming  and  need  a  lot  of  storage 
space  [Grin'iths  o']  but  they  can  accelerated  by  using  in  addition  an  estimation  of  the  intersection 
region  which  yields  to  the  so  called  ' dhide- and- conquer'  algorithms.  For  B-spLine-Bezier  repres¬ 
entations  the  estimation  of  those  parts  that  do  not  participate  in  the  intersection  can  be  done  by 
usmg  the  convex  hull  property  [Lane  et  al.  '50].  [Peng  '8-f]  (see  also  part  1  of  this  paper)  or.  rougher 
but  more  easily  and  quickly  handled  by  min-max  boxes  (see  part  II  of  this  paper).  .-Vn  estunauon 
using  min-max  bo.xes  can  also  be  done  for  non-B-spline- Bezier  representations  [Koparkar.Mudur 
■55], 

.A  disadvantage  of  the  subdivision  and  even  of  the  more  advanced  divide-and-conquer  algorithms 
again?t  the  algebraic  based  intersection  algonthms  migkt  be  that  they  are  more  time-consummg 
^^ecause  of  the. subdivision  process  [Sederberg, Parry  '86],  But  the  great  advantages  of  the  divide- 
and-conquer  algorithms  are  that 

•  fhey  are  very  user-fnendly  -  no  worry  about  'suitable'  starting  points, 

•  they  find  independently  -  that  means  without  any  interactive  disruption  to  the  user  -  all  inter¬ 
section  points  w  ithin  the  specified  tolerance, 

•  they  can  be  formulated  easily  for  arbitrary  polynomial  degree  and  for  non-rational  and  rational 
representations,  and 

•  they  are  numerically  very  stable  because  of  the  extraordinary  numerical  properties  of  the  Bernstein 
polynomials  [Farouki.Rajan  '57],  [Sederberg.Pany  '86], 


Because  of  these  favorable  properties  of  the  intersection  algorithms  based  on  the  B-spline-Bezier 
representation  using  divide-and-conquer  methods  and  because  of  the  reason  named  above  we  w  ould 
like  to  have  also  a  self-intersection  algorithm  of  this  kind.  The  existing  curve-curve  algorithms  can 
not  be  used  directly  by  doing  the  curve  input  twice,  because  in  this  case  the  divide-and-conquer 
method  will  fail  in  the  sense  that  no  elimination  of  curve  parts  that  do  not  participate  m  the  self- 
intersection  will  be  possible.  Furthermore  the  fmal  calculation  of  the  self-intersection  points,  done 
by  intersecting  straight  line  segments  defmed  by  the  control  polygon  will  also  fail  by  doing  the  same 
polygon  mput  twice. 


The  only  self-intersection  algorithm  for  B-spline-Bezier  representations  I  know  was  given  in 
[Tiller, Hansen  '84],  They  calculate  self-intersections  of  (rational)  B-spline  cuives  in  a  two  step 
procedure.  First  they  fmd  the  intersections  of  the  control  polygon  with  itself  and  then  they  use  an 
iterative  method  (e.g.  Newton)  to  improve  the  approximate  solution  found  in  step  one.  They  know 
that  this  method  can  fail,  because  a  cuive  can  have  a  loop  even  though  its  control  polygon  has  no 
sell'-intersection,  but  by  using  control  polygons  which  approximate  their  curves  very  closely,  i.e. 
building  up  the  curve  by  a  'large'  number  of  segments,  they  try  to  make  sure  to  be  'on  the  safe  side'. 
.Although  Tiller  and  Flansen  are  working  with  B-spline  techniques,  their  algorithm  do  sen  t  belong 
to  the  powerful  class  of  the  divide-and-conquer  algorithms  because  their  algorithm  dosen't  use  the 
typical  kind  of  strategy  of  the  divide-and-conquer  algorithms  for  the  evaluation  of  the  self- 
intersections. 


The  algorithm  presented  here  is  a  user-friendly  divide-and-conquer  algonthm  for  finding  all  the 
self-intersection  points,  including  their  parameter  values,  of  a  parameterized  non-rational  or  rational 
curve  of  arbitrary  degree  in  Bezier  representation.  For  the  creatmg  of  the  algorithm  the  geometnc 
relationship  between  the  curve  and  its  defining  control  polygon  was  fully  taken  into  account.  By 
alternately  subdividing  the  Bezier  polygon  and  esiimatmg  the  self-intersection  regions  the  self- 
mtersection  points  are  fmally  approximated  by  straight  Une  intersections  of  the  refmed  Bezier 
polygon.  In  addition  to  the  convex  hull  property  and  the  approximation  property  of  the  Bezier 
polygon  the  algonthm  is  based  on  a  very  mtuitive  angle  critenon  which  is  together  with  the  convex 
hull  property  used  for  estimatmg  the  selJ'-intersection  region  of  the  curve. 

Because  a  cune-curve  intersection  algonthm  is  an  important  part  of  the  self-intersection  algorithm 
o!  part  in  of  the  paper,  and  because  the  final  calculation  of  the  self-intersection  points  and  its  pa¬ 
rameter  values  is  done  in  the  same  way  as  m  the  curse-curve  algonthm,  a  short  explanation  of  a 


divide-and-conquer  algorithm  for  calculating  the  intersection  points  of  two  parameterized  non- 
rational  or  rational  curves  of  arbitrarv’  degree  in  Bezier  representation  is  given  in  part  II,  The 
cune-cune  algorithm  described  there  differs  from  the  classical'  one  mtroduced  by  Lane  [_Lane  a 
al.  fd]  in  some  details' ,  mainly  in  the  concept  of  the  'control  unit'  and  in  the  fmal  calculation 
of  the  intersection  points  and  its  parameter  values. 

Part  I\'  finally  mcludes  a  short  description  of  how  to  calculate  the  self-intersections  of  a  Bezier 
spline  curve. 

.-Ul  algorithms  are  written  for  planar  curves,  but  for  the  e.xtension  to  spatial  curves  only  'a  third 
eiiuation  for  the  z-coordinate  '  has  to  be  added  everywhere  where  coordinates  have  to  be  evaluated. 
The  paper  starts  with  some  introductory  words  on  the  Bezier  representation  of  (planar)  curves. 


/.  Bezier  Curves 

A  (planar)  Bezier  curve  is  defined  by 

m 

B(u)  =  y 

fc-O 

where  b*  =  (ar*  ,>'J  e  RL  ue  [0,1]  and 

5'"/  \  I tn\  k,,  ,m—k 

ki^^)  = 

are  the  (ordinary)  Bernstein  polynomials  of  degree  rn  'inu.  The  coefficients  b»  e  are  called  Bezier 
points.  They  form  in  their  natural  ordering  given  by  their  subscripts  the  vertices  of  the  so  called 
Bezier  polygon  (see  Figure  1). 

It  is  possible  to  build  up  complex  Bezier  spline  curves  from  a  number  of  Bezier  curve  segments. 
The  conditions  for  C  continuity  of  adjacent  curve  segments  can  be  found  m  [Boehm  et  al  (?•/], 

The  Bezier  description  of  a  curve  is  a  very  powerful  tool  because  the  expansion  in  terms  of 
Bemstem  poivnomials  yield  to  a  geometric  relationship  between  the  curve  and  its  defining  Bezier 
points.  For  example: 

•  the  Bezier  polygon  gives  a  rough  impression  of  the  Bezier  curve  (sec  Figure  1), 

•  the  cur.e  has  its  endpomts  at  b,,  and  b,  with  tangent  vectors  defmed  by  bg,  b,  and  by  b,,  b,., 

1  see  Figure  I ). 

•  consex  'hull  property:  the  Bezier  curve  lies  completely  within  the  convex  hull  of  its  Bezier  polygon 
('ce  Figure  2). 

•  the  cur.  e  pouit  Bfuj).  for  any  e  [0.1]  can  be  computed  by  repeated  de  Casteljau  steps  by  the 
recursion  formula 

bf(i^)  =  (1  -  uo)b^~'(uo)  +  Uob^+,(Uo) 
where  b;  s  b.  and  B(i^)  =  by  (see  Figure  3). 

The  point  n  =  u,  subdivides  a  Bezier  curve  uito  two  C”  continuous  segments.  Each  segment  is  again 
a  Bezier  cur.e  of  the  same  degree  as  the  ongmai  one.  The  Bezier  pomts  of  these  two  segments  are 
b'.produc'.i  of  the  de  Casteljau  construction  for  the  evaluation  of  the  point  Bm).  They  are  given 
by  bf  and  b;  (4  =  0,  .,  .m).  The  subdivision  process  may  be  repeated  yielding  a  sequence  of 
polygons.  For  this  sequence  of  polygons  we  have  the  important 

•  approximation  property:  if  the  u,  are  dense  m  [0.1]  the  sequence  of  polygons  converges  to  the 
curve. 

Figure  4  illustrates  how  the  curve  can  be  fi.xed  usmg  the  appro .ximation  and  the  convex  hull  prop- 


.  IT*  ’ J<  ~Jf. 


Figure  4.  fixing  the  curve  by  the  approximation  and  the  convex  hull  property 


A  raiionsU  (planar)  Bezier  curve  can  be  defined  by 

R(«)  =  V  bfciR^(u) 

r  asii 

'.'.here  bi  =  (.t»  .y*)  e  ue  [0.1]  and 


R'iu) 


/?■■  B"(u) 


are  the  rational  Bernstein  polynomial  of  degree  rr;  m  u  with  weights  e  R  [/’teg/ 

f’lCJre  5  compares  the  (ordinary)  Demsiem  polynomials  B.'iu)  and  the  rational  Bernstein 

polynomials  RTiu)  with  /I,  >  U  for  all /c- 

If  ••'.e  demand  /J*  >  0  for  all  k  we  have  aU  the  properties  and  algonthms  for  rational  Bezier  curves 
■Ahich  we  have  for  orduury  i.e.  non-rational  cur-es  [farm  .sj],  [_Tiikr  SJ],  therefore  there  is  ui  this 
cj'C  no  principle  difference  between  a  eurve-eune  resp.  a  curve  seU'-interseetion  algonthrn  for 
non-rational  and  for  rational  Bezier  curves. 


Figure  S.  ordinary  and  rational  Bernstein  polynomial  of  degree  four,  (  . =  (1.  3.  2.  5.  I  ) 


//.  Curve-Curve  Algorithm 

ITie  underlying  idea  of  the  cune-cune  algorithm  is  to  deal  with  the  Bezier  polygon  instead  of  the 
eur.e  desenption  itself,  using  the  relations  between  polygon  and  curi'e  mentioned  above. 

I  he  program  of  the  algonthm  is  to  subdivide  both  cur\'es  repeatedly  which  yields  at  the  same  time 
to  a  subdivision  and  refmement  of  the  polygons.  This  is  done  until  a  fine  polygon  structure  is  ob¬ 
tained  and  the  curves  can  be  appro.vimated  well  by  the  polygons  defined  by  these  subdivisions.  Thus 
procedure  reduces  the  problem  to  a  number  of  straight  line  intersections  that  can  be  handled  easdy. 
Because  subdividing  the  whole  curves  in  each  algorithm  step  is  relatively  time-consuming  and  needs 
a  lot  of  storage  space  m  addition  an  estimation  of  the  mtersection  region  is  done. 

The  algonthm  consists  of  four  mam  parts  (Figure  6),  they  are  described  now. 

•  First,  the  intersection  area  is  estimated.  Using  a  coarse  but  very  quick  estimate  of  the  possible 
intersection  regions  of  the  two  curves  those  parts  of  the  curves  that  do  not  participate  m  the 
intersection  will  be  elimmated  as  early  as  possible  in  the  algontlim. 

•  Second,  refmement  occurs  by  subdividing  the  Bezier  polygons.  Except  at  the  beginmng.  the  al- 
contlim  subdivides  not  the  whole  Bezier  polygons,  but  only  those  pans  whose  corresponding 
cur.  e  ponions  might  participate  in  the  mtersection.  -An  adaptive  subdivision  is  done  to  detect  the 
separation  of  regions  of  the  two  curves  that  do  not  intersect  readily. 

•  T  turd,  the  intersection  pomts  are  calculated  by  mtersecting  the  Bezier  polygons  of  the  curve 
sub-egments  of  possible  mtersection.  Part  three  also  calculates  the  parameter  values  of  the 
;nter-ection  points. 

•  Tounh.  error  values  are  calculated,  tolerances  are  checked,  this  part  of  the  algorithm  is  the  con¬ 
trolling  umt  of  the  algonthm  and  is  very  important  tor  deaimg  with  difficult  and  complicate  cases. 

Beside  drawmg  parameters  for  creating  the  plot  output,  the  input  of  the  algonthm  consists  of  the 
polynomial  degrees  iM  and  m)  and  of  the  Bezier  pomts  of  the  two  Bezier  curves  (B(7T  and  bini. 
furthennore  of  an  error  tolerance  value  to  determme  the  accuracy  needed.  Pre-settmgs  for  con¬ 
trolling  the  algonthm  can  be  specified  m  the  program  too. 

The  first  step  of  the  algonthm  is  to  subdivide  the  two  curves  simultaneously  forming  two  new 
sub'eginents  on  each  curve.  A  nun-max  box  defmed  by  the  maximum  and  minimum  -v  and  _v  co¬ 
ordinates  of  the  curve  segments  defmmg  Bezier  points  is  built  for  each  segment.  The  boxes  of  the 
f.vo  evrves  are  then  compared  with  each  other  (a  companson  usmg  min-max  boxes  instead  ol  the 
convex  hulls  is  rougher,  but  much  more  easily  handled  and  quickly  practised).  Those  subsegments 


6 


whose  boxes  do  not  intersect  any  box  of  the  other  cune  will  no  longer  be  considered.  Only  those 
subsegments  whose  boxes  can  not  be  separated  from  that  of  theu  avals  will  be  dealt  with  funher 
(Figure  ').  For  this,  Bezier  points  of  pairs  of  interfering  subsegments  of  different  curves  will  be 
pres  ided  ssith  an  subsenpt.  called  interference  index'.  By  this  a  list  of  pairs  of  segments  of  duferent 
curves  uhich  might  interfere  is  created.  In  the  following,  Bezier  points,  i.e.  segments  of  the  'ume 
interference  index,  vsiU  all  go  through  the  algorithm  subroutines. 

The  de  ("asteljau  subdivision  process,  the  mm-max  box  formation  and  the  separability  test  are 
connected  by  an  algonthm  loop,  which  will  be  done  as  often  as  is  required  by  the  lev  el  of  accuracy 
needed.  .-Viter  each  subdivision,  two  new  subsegments  are  formed,  each  correspondme  to  a  smaller 
convex  hull,  ^^■flen  more  and  more  subdivisions  are  done  each  convex  hull  becomes  smaller  and 
smaiier.  while  the  cur.e  topology  near  the  intersection  is  reasonably  closely  approximated  b>  the 


Figure  7.  estimating  the  intersection  region  using  min-max  boxes 

•Ml  subsegments  which  might  participate  in  the  intersection  go  through  the  third  part  of  the  algo- 
r.thm:  the  section  that  computes  the  intersection  points  and  the  parameter  values  of  the  intersection 
points  what  is  be  done  in  the  follow iniz  wav. 


Let  a  subsegment  of  the  first  cup. e  B(7^  of  degree  .Vf  and 


B;  =  [BXj.BVj) 


/  =  0 . \f 


be  the  Bezier  points  of  B,(t)  and  let  a  subsegment  of  the  second  curve  b(t)  of  degree  m 

and 


b/e  =  ibx^  .  bvft) 


k  —  0 . m 


be  the  Bezier  points  of  b,^(T). 

The  polygon  legs  defined  by  the  Bezier  points  are  given  by 


Gj  =  By  +  T,Sj 


y  =  0 . V/  -  1 


where  G,  =  (f/A’. ,  G  S.  =  (SA' .  S  K,).  S.  =  B..,  B.  and  7  e  [0,1]  and  similar  for  g.. 
If  G  and  g*  intersect  in  P  (Figure  S)  i.e. 

G(fy=f(P))  =  P  =  g*{t*  =  4(P)) 


we  have  for  the  parameter  values 

f,(P)  = 


svi^iBX,  -  bxfj  -  sx/^iBYj  -  bvj^) 


4lP)  = 


SYj(BXj  -  bxi,)  -  SXjiBYj  -  by^) 


where 


■\(t  “  ^-^It  Vit 


7.  P)  resp.  h(P)  are  parameter  values  with  respect  to  the  polygon  legs  G,  resp.  g.  but  because 
the  de  Casteljau  refmement  is  always  done  for  0.5  we  also  know  the  parameter  value  7-  of 
B,  =  B,.(0)  and  the  parameter  value  of  b,  =  b;f(0)  so  that  the  parameter  values  7(P)  and  nPi 
ot  the  intersection  point  P  with  respect  to  the  parameter  intervals  of  the  onginally  given  Bezier 
car.es  can  be  calculated  by  (Figure  7) 


7(P)  =  7, 


7,  +  (7^^,  -  7.)  7.(P) 


and  suTular  for  .'(Pi,  where  i  is  the  number  of  subdivisions  and  7  are  the  parameter  values  given 
to  the  Bezier  points  B  of  the  Bezier  polygon  of  BAt).  The  7,  (and  so  the  t,  given  to  the  bj 
can  be  defined  m  dilTerent  ways  for  example 


by  an  equidistant  measure 


r  = 


by  an  chord  length  measure 

1 

1=IJ 

by  an  geometric  average  measure  of  7'  and  7 


w/iere  L  = 


M-\ 

_  V 


B  -  B 


8 


Figure  8.  calculating  of  parameter  values  of  the  intersection  points 


As  a  measure  of  error  w  e  can  use  the  distances 

Rsb  =  llB(r(P))  -  b(r(P))|| 


Rsp  =  llB(r(P))  -  P|| 
R^P  =  i|b(r(P))  -  Pll 


Per  default  a  minimum  number  of  de  Casteljau  subdivisions  will  be  done  before  part  three  will  be 
started  (loop  1).  If  the  accuracy  needed  is  0.002  for  example  the  pre-setting  has  to  be  6  (see  table 
1)  an  this  will  yield  in  almost  every  example  to  an  accuracy  of  about  0.002,  if  in  some  complicated 
case  not,  the  control  unit  will  effect  to  do  as  many  additional  subdivisions  as  needed  for  the  speci¬ 
fied  accuracy  (loop  2). 

When  the  two  curves  intersect  in  a  very  small  angle  or  do  not  intersect,  but  come  very  close  together 
part  three  might  calculate  more  intersection  pomts  as  two  curves  of  degree  M  and  m  can  produce 
or  might  calculate  ( pseudo- lintersection  points  lying  very  close  together  in  parameter  space  which 
has  to  be  checked  (the  statement  of  the  parameter  space  criterion  is  stronger  than  an  statement  of 
an  analog  coordinate  space  cnterion).  In  both  cases  the  control  unit  will  also  effect  to  do  as  many 
additional  subdisasions  as  needed  for  clarifying  the  situation. 

The  repeatedly  done  polygon  reffnement  iiutialised  by  these  criterions  will  be  stopped  in  different 
wa>  s:  lirst.  d  the  result  has  the  accuracy  needed,  second,  there  is  a  default  of  an  upper  boundary  for 
the  number  of  de  Casteljau  subdivisions  and  third,  there  is  a  default  of  an  maximal  (possible)  ac¬ 
curacy  This  default  value  is  dependent  on  the  initialization  of  the  vanables,  e.g.  real  or  double 
precision  real  and  of  the  machine  accuracy  for  each  kind  of  Initialization. 

Fmall>  the  control  unit  checks  if  the  distance  between  intersection  pomts  in  coordinate  space  is  less 
than  a  specilied  tolerance.  If  yes.  an  mtersection  pomt  is  defined  by  the  arithmetic  average  of  these 
points 


Examples 

Table  I  lists  the  ma.ximal  error 

R  =  tnax  (  Rg^  .  Rgp  .  R^p)  0 

as  it  depend  upon  an  mcreasmg  subdivision  factor  for  the  examples  1  to  6  for  equidistant 
parameterization  lor  wluch  we  got  the  best  results. 


parameter  values  and  x-j -coordinates  of  the  intersection  points 


PX 

PY 

r(P) 

t{P) 

-3.^4353 

1.49822 

0.23120 

0.27305 

-2.9:393 

1.500S6 

0.29330 

0.32148 

■0.7"3:5 

1.49989 

0.44827 

0.45409 

0."'325 

1 .49989 

0.55173 

0.54591 

2.92393 

1.511086 

0.70670 

0.67852 

3.64353 

1.49822 

0.76880 

0.72695 

III.  Self-Intersectii 

It  is  not  possible  to  calculate  t 
of  part  II  by  doing  the  cur.  e  in| 
al'.'.  a;. '  be  poMtne  'O  that 
intors.vtion  is  possible.  Furthci 
aJuitK.'iiai  cntenon  is  necessary 
^^'hat  'ae  ould  like  to  lu'.e  is 
ueilrung  Beeier  pouits  i  e.  its  B 
the  cor.'.ew  huh  prcperts  .  This 
'ifaation  IS  compheated  by  the 

•  It  IS  possible  that  the  Bezier 
intersection  isee  Figure  9) 


but  on  the  other  side  even 


no  self-intersection 


Figure  13.  Bezier  curve  with  self-intersection 


“2 


Figure  15.  j*  with  different  orientation 


Because  of  the  de  Casteljau  construction  u  hich  creates  in  evers  step  a  convex  combination  of  the 
and  because  of  the  appro.ximation  property  mentioned  in  part  I.  the  sum  21  ixj  is  equal  to  the 
sum  of  the  amount  of  the  rotation  an  de  of  the  tangent  vector  of  the  Bezier  curve  if  the  orientation 
of  the  rotation  angles  of  the  Bezier  polygon  legs  is  the  same  in  evens'  inner  Bezier  point  (see  Ficure 
I4i.  But  the  sum  of  the  amount  of  the  rotation  angle  of  the  tangent  sector  of  the  Bezier  cur. e  is 
smaKer  than  21  ^  'he  orientation  of  the  rotation  angles  of  the  Bezier  poKgon  legs  is  not  the 

sam.e  in  eser.  inner  Bezier  point  (see  Figure  15'.  that  foUosvs  from  the  smoothing  properts  of  the 
de  Casteljau  jubdisision  process  together  with  the  approximation  property  mentioned  m  part  I, 

So  sse  base  the  statement  that 


■  the  sum  of  the  am.ount  of  the  rotation  angle  of  the  tangent  vector  of  the  Bezier  curve  is  alssass 
smaller  or  equal  the  sum  2i  '  ol  the  amounts  of  the  rotation  angles  of  the  Bezier  pols  gon  legs. 

By  combining  the  two  ■  statements  sve  get  the 


•  ansle  criterion:  The  sum  21  of  the  amounts  of  the  rotation  angles  of  the  Bezier  polygon  legs 
is  greater  than  t  if  the  Bezier  curve  has  a  self-mtersection. 


Tor  the  algorithm  sse  ss  Ul  use  the  contraposition  of  the  criterion. 


•  on'ile  criterion:  A  Bezier  curve  has  no  self-intersection  if  the  sum  2i  l®*i  of  the  amounts  of  the 
rotation  angles  of  the  Bezier  polygon  legs  is  smaller  or  equal  than  n. 

By  this  ss  e  have  a  s  ery  simple  geometric  criterion  for  deciding  svhether  a  curve  has  a  self-intersection 
or  not  and  for  the  elimination  of  curve  parts  that  do  not  participate  m  a  self-intersection.  VVhat  sie 
base  to  do  is  to  calculate  the  sum  21 1**1  of  ^he  polygon  angles  and  compare  with  t  .  If  s\e 
base  2i  ■*.!  ^  knosv  that  there  is  no  self-intersection  of  the  curve  (Figure  16.1  and  16,2).  but 
if  21  ®.i  >  ^  the  curse  might  have  a  self-intersection  (Figure  16.3  and  16.4).  For  clarify  mg  sse 
subdiside  usmg  de  Casteljau  and  check  the  smaller  parts  again  against  the  angle  critenon. 


Figure  16.  the  angle  criterion 
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To  buiJJ  up  a  SL'lf-intcrsection  algonthm  the  idea  of  the  angJe  test  has  to  be  combined  u  uh  the  idea 
of  the  rrun-ma.x  bo.\  test.  I  his  is  done  in  the  foUoisme  way. 

The  aleonthm  consists  aeaui  of  the  four  main  parts  of  Figure  6.  But  part  one  of  the  aigonihm  for 
estumatmg  the  self-intersection  region  of  the  curse  consists  now  of  tsso  diilereni  tests,  tlie  mm-max 
bo.x  test  -md  the  angle  test.  1  igure  In  g:\es  the  e.xample  of  an  subdis'ided  Bener  cur.e  haiine  -e- 
\erai  seit-mtersections.  As  '.se  can  see,  there  is  a  subsegment  (subse'gment  B<)  with  ^elf-mtersecuon 
(point  P  i,  th.re  are  two  subsegments  with  common  boundary  point  (subsegments  B,  and  B,)  cre¬ 
ating  the  'clt-mtersecticn  point  P;  and  there  are  subsegments  (subsegments  B|  and  B, i  which  are 
not  connected  to  each  otlier  but  create  the  self-intersection  point  P,  of  the  Bezier  cur.e  Bi  T'.  I  o 
di'tir.guish  between  these  three  ddferent  cases  and  for  controlling  the  algorithm  we  mtroduce  a  so 
called  zfn.r  .'.■tttt’.v 

1  he  -ell-intersection  of  a  segment  of  genus  one  that  means  a  segment  of  case  one  has  to  be  checked 
b>  usmg  the  angle  criterion.  If  the  angle  test  is  positive  e.g.  ^  in,;  >  it  a  refinement  has  to  be  done 
to  cLnfy  the  situation.  The  retlnement  of  a  genus  one  segment  produces  two  subsegments  of  genus 
one  and  one  pair  of  subsegments  of  genus  two. 

A  parr  of  subsegments  of  genus  two  that  means  subsegments  with  a  common  boundary  pomt  hate 
also  to  be  checked  agamst  the  angle  cnterion  but  now  the  angle  sum  of  both  polsgons  has  to  be 
cA'ulated.  The  min-ma.x  bo.x  entenon  can  not  be  used  because  of  the  common  boundarv-  pomt 
of  the  two  segments.  If  the  angle  test  is  positive  a  refinement  has  to  be  done  for  both  segments,  it 
produces  one  pair  of  subsegments  of  genus  two  and  three  pairs  of  subsegments  of  genus  tliree. 
Subsegments  of  genus  tliree  can  be  dealt  with  as  in  the  curve-curve  algonthm  ot  part  11  i.e.  for 
calculatmg  the  sell-mterscction  pomt  P,  of  Figure  16  we  do  need  only  the  min-max  box  test  not 
the  angle  test  because  the  retlnement  of  pairs  of  segments  of  genus  three  can  produce  pairs  of  sub- 
segments  of  genus  three  only  and  no  (pairs  of)  subsegments  of  genus  one  or  two. 


Figure  17.  possible  subsegment  configurations  contributing  to  the  self  intersection 

In  the  flr^t  step  the  algonthm  has  to  deal  only  with  one  curve  segment  of  genus  one  -  the  Bezier 
cur.  e  'Cgment  which  has  to  be  checked  for  self-intersections.  If  the  angle  test  is  positive  a  refmement 
has  to  be  done,  so  that  m  the  second  step  the  algonthm  has  to  deal  with  two  subsegments  of  genus 
one  and  one  pair  of  subsegments  of  genus  two  .ind  the  result  of  this  step  might  be  subsegments  of 
genus  one  and  pairs  of  subsegments  of  genus  two  or  three.  SVlien  more  and  more  subdivisions  are 
done  not  only  each  convex  hull  becomes  smaller  and  smaller  but  because  of  the  appro.ximation 
propery  of  the  Bezier  polygon  also  the  ancle  -ume  of  each  subsegment  becomes  smaller  and 
smaller  lo  that  after  an  initial  increase  of  (pairs  ot'i  'ubsegments  of  genus  one  and  two  the  number 
of  (pairs  of)  subsegments  of  genus  one  and  two  decreases  very'  fast  until  there  are  only  pairs  of 
subse'gments  of  genus  three.  From  thus  m.omcn:  on  the  self-mtersection  algonthm  works  m  the 
same  way  as  the  curve-curve  algonthm  dc'cnbcd  m  part  II  of  the  paper.  That  also  means  that  p.irt 
two  and  part  three  of  the  algonthm  -  the  subdr.  mon  of  the  curve  in  the  aim  of  refinement  and  the 
calculation  of  the  self-intersection  pomts  and  parameter  values  -  is  done  in  exactly  the  same  way 
as  described  m  part  II. 


:o 


The  control  umt  also  works  as  m  part  II  Jescabted  except  that  it  checks  in  addition  d  the  subdi- 
\ided  control  pohgon  turns  through  ISO  degrees  at  a  subdivision  point  which  implies  a  cusp  at  tius 
pomt  (Example  4). 


Examples 

Table  2  lists  the  maximal  error  given  by  >g)  part  II  as  it  depend  upon  an  increasmg  subdivision 
iactor  ior  the  examples  1  to  12  for  equidistant  parameterization. 


subdiMsion  factor 

4 

5 

6 

7 

8 

Example  1 

0.2S064 

0.01093 

0.00568 

0.00269 

O.iJiJl  10 

Example  2 

0.08635 

0.03628 

0.01297 

0.00189 

0.uiHj"5 

Example  3 

0.20799 

0.09101 

0.02185 

0.00619 

0.00 101 

Example  4 

Example  5 

0.02390 

0.00445 

0.00156 

0.00007 

0.i.)|.i(i(j3 

Example  6 

0  03326 

0  01207 

0.00236 

0.00080 

O.OiJO'iO 

E.xample  ~ 

0.04659 

0.01579 

0.00055 

0.00027 

0  00012 

Example  8 

0.06235 

0.02459 

0.00620 

0.00150 

0.00040 

Example  9 

0.10368 

0.01132 

0.00434 

Example  10 

0.10184 

0.01123 

0.00427 

0.00159 

0.00032 

E.xam.ple  11 

E.xample  12 

0.07962 

0.00794 

0.00332 

0.00124 

0.00025 

Table  2.  R  for  equidistant  parameterization 


Because  of  the  bad  character'  of  the  two  cusps  appearing  in  Example  II,  this  example  requires 
more  than  S  subdivisions  for  the  decision  if  the  curve  has  self-uitersections  or  cusps. 


21 


(Ji’liiyDO 

5  1 

parameter  'alues  and  x-> -coordinates  of  the  intersection  point 


parameter  values'  and  x-y-coordinates  of  the  intersection  points 


Example  6 


Example  8 


-Q 


O 


PX 

PY 

t,(P) 

1.5'-')322 

|)4185‘) 

parameter  values  and  x-y-coordinates  of  the  intersection  point 


'  X 

.!■: 

0  ij 

1  .  , 

0  0 

M  .  1 

0  0 

'  1  1 

14  0 

\  '  1 

14.0 

0  0 

^  0 

Be/ier  points  of  b(/) 


•0,65030 

he 

"3 

i  1017730 

i),5S924 

0,00000 

0  22167 

0.77S33 

IV.  Self-Intersections  of  Spline  Curves 

^■o^ma!;>  are  not  reaL;^  interested  in  <inde  [kvier  cur/e  se^nents  but  m  B-splme  resp.  Bev.ier 
'piine  wi.r. es  oon'istine  v't  'O'. eral  our. e  segments.  Because  a  B-splme  curve  can  be  redefined  m  .1 
Beeier  ‘onn  b;.  u-une  ti’.e  (  Klo  aijtinthm  addine  multiple  knots  m  one  pass  \_Conen  et  al.  >"]  or 
isy  U'lne  tne  ^orr.putation.iay  more  eilicient  Boehm  ateontiim  adding  tlte  multiple  knots  one  b\ 
one  S  '.  >_'j,  'Oit-ir.tersections  ot  B-spline  and  01  Bezier  spline  cur.es  can  he  calculated 

using  ti'.e  aigor.tr.tr,'  ol'  part  II  and  111. 

The  segments  ai  ot  the  Bezier  representation  ot  the  splme  cur.e  might  be  snven  bv 

V.  here 


a  =  {I  —  u)  u  ,  0  <  u  <  I  ,  A'  =  0 . 

i  e.  the  splme  curve  is  defmed  with  respect  to  a  partition  of  the  domam  space  by  knots' 

Zq  <  Z|  <  ... 

The  self-mtersection  points  ot  a  spline  curve  can  be  calculated  by  domg  the  curve-curve  mtersecticn 
algonthm  for  all  pairs  of  segments  and  with  K  *  K  and  by  domg  the  curv  e  self- 
mtersecticn  algonthm  for  all  segments  .  While  the  aigonthms  of  part  11  and  111  calculate  the 
parameter  val  les  of  the  seLf-mtersection  pomts  with  respect  to  the  local  coordmate  domam  [0.1] 
we  also  Lnow  •  because  of  the  Imear  relation  between  z  and  u  •  the  z  parameter  values  of  the 
seil'-mtersection  pomts. 


Remark 

llus  study  was  done  as  a  pre-study  for  the  creaimg  of  a  surface  self-mtersection  algonthm  for 
pz.rami-etcrLzed  surt'aces  m  Bezier  representation.  The  surface  algonthm  is  descnbed  m  the  paper 
Self- Intersections  of  Parametric  Surfaces  .  TecTmical  Report  »  \PS-53-88-002.  Naval  Postgraduate 

School.  Monterev  (lOSSi. 
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